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Ordinary Bessel Functions (DOI OR SADOl) 

Closed with program parameter 

257 

0, 1, 2, 3; 6, 7, 8* 9 

About k decimal places (see below) 

At most 100 ms. (see below) 

-M 
Enter with x -2 in A, M > 



p 


XY MF 
50 pF 


p+1 


26 gF 



(this routine at g) 



X * 5 -» Functions of first kind (J:-) only will be computed 

X = J -4 Functions of both first and second kinds (Yn often written N ) 

n 

will "be computed 
1 = 1 -4 Functions of zero order §T *, Y ) only will be computed 
Y * 2 <-» Functions of first order (JT* , t-A only will he computed 
y = -» Functions of both zero and first orders will toe computed 

More explicitly: 



XY 
5k 
52 
50 
Jk 
J2 
JO 



XnjB. computed 

O 

J (x>* J x (x) 
J (x>, Y Q (x) 
J x (x), Y x (x) 

j q (x), ^(x), y 6 (x), Y x (x) 



19 /- 

The results are scaled by 2 and are left in locations o - 9 as follows 

J (x) .2 9 in loo 6 
o 10 

Jjfr) .2" to lac 7 



- 2 - 



Y (x) .2~ 19 in loc 8 
Y 1 (x) .2~ y in loc 9 



Illiac will stop if x <C 0, and upon starting again will transfer control 
to the right hand order at the location specified "by the preset parameter 
S3- 

NOTE 



,-M. 



The integer M need not be the smallest integer such that X.2 < 1. 

DESCRIPTION If x ^16, the following asymptotic expressions are used: 

2 r . . /. . *% . 1 



J e (x) x \j— [sin (x + £) + jg sin (x - |)] 
J x (x) £ \;f£ t sin ( x - f) + Hx* Sin ^ x + ¥^ 
Y o (x) - ^ Csin (x ' ^ * S Sin (X + £ )] 

No natter what entry is used;, if x ^ 16, then all four functions are computed* 
The time is about 25 ms . At x - 16, where these f orsaulae give the worst 
results, the error is less than .0001. 
If < x < 1.6 , the following expansions are used 

N /x s2n 

J (x) = Z n (-l) n { 2 } 

o x n=0 v *- 

(ni f 

N ( X ) 2n 

J l (x) = f n?0 ( - 1)n nl (n+l)l 



o N 

Y o (x)=|[^ + 0,|) Jo W-nSl 



/X v 2n 
(-if ( 2 ) 



(nl)' 



— (1 + 1+ •••+!)] 



V x) = I [(tf+Mf) J i (x) - 1 ] 



- 7 to (-D° ( S (i + | + • 

nl (n+1) 1 



1 _ 1 
+ n + 1 + 2 + 



n + 1 



)] 



3 - 



where N is large enough in each ease so that the last term is less than 2 
The time to compute J (x) or J. (x) is at most kO ms . depending on the size 

of x (less time for smaller x)» The time to compute J (x) and Y (x) or 

° -16 ° 

The error is less than 2 



-16 



J. (x ) and Y, (x ) is at most 100 



ins 



If J (x) or Y^Cx) with n 4 s 0, 1 are desired j, use the recurrence relations: 

n n ■ 

2n 



T- J n W = J n + 1 W + Vl< x > 

V x > - W x > + W x > 



2£ 
x 



For x < Q° 9 



J (0) * 1 

o v ' 

J x (0) = 



Y Q (0) = - ^ 

Y 1 (0) m -^ 

We also have 

J (- X ) = J ( x ) 


J x (-x) = - J x (x) 

Y (-x) = Y q (x) + (2n+l) *i J q (x) 

Y-^-x) = Y x (x) - (2n+l) *i J^x) 

The branch point at the origin produces the multiple values of 

Y (-x) and Y- (-x) 
o l v 

Library routines SI (Natural Logarithm) and T5 (Sine - Cosine) are included 

in this routine and may "be used by other parts of the program . SI starts at 

word g + 191 and T5 starts at word g + 236, where g is the location of the 

first word of this routine* 
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